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LONG  TERM  VARIABILITY  OF  WINDSPEED 


INTRODUCTION 

Measurements  of  ambient  ocean  acoustic  noise  at  frequencies  from  a  few  hundred 
hertz  to  several  kilohertz  generally  show  strong  correlation  with  wind  speed.  Ambient  noise 
is  correlated  with  surface  wave  height,  also,  but  to  a  lesser  degree,  and  thus  wind  speed 
appears  to  be  the  single  most  important  weather  variable  needed  in  making  predictions  of 
ambient  noise.  In  support  of  the  need  of  the  environmental  acoustics  community  to  be  able 
to  predict  ambient  noise  variability,  this  report  describes  the  statistical  variability  and  fluc¬ 
tuation  spectrum  of  long  term  records  of  wind  speed,  obtained  from  two  weather  stations. 
A  brief  review  is  given  also  of  results  from  previous  studies  of  wind  speed  fluctuations.  The 
present  results  are  used  in  part  to  extend  the  frequency  range  of  previously  reported  wind 
speed  fluctuation  spectra.  Time  scales  of  interest  here  are  those  corresponding  to  synoptic 
variability,  which  has  periods  of  days  to  months,  and  to  seasonal  or  yearly  variation.  Brief 
attention  is  paid  also  to  records  of  sea  level  atmospheric  pressure.  This  variable  shows  cor¬ 
relation  with  wind  speed  and  may  be  of  some  use  in  noise  prediction. 

The  data  that  form  the  basis  of  this  report  include  values  of  surface  wind  speed  re¬ 
corded  every  12  h  over  the  2-yr  interval  from  1  June  1974  to  31  May  1976  at  Bear  Island 
(World  Meteorological  Organization,  WMO,  Station  1028,  74°N,  19°E)  and  at  Jan  Mayen 
(WMO  Station  1001,  71°N,  8°W).  Bear  Island  is  located  in  the  Barents  Sea,  and  Jan  Mayen, 
in  the  Norwegian  Sea.  A  second  and  smaller  data  set,  consisting  of  wind  speed  values  every 
3  h  from  1  July  1974  to  31  January  1975  at  Bear  Island,  was  also  available  and  was  used  to 
extend  the  high-frequency  end  of  the  fluctuation  spectrum.  These  two  stations  were  chosen 
in  order  that  the  results  could  add  to  other  Naval  Research  Laboratory  studies  of  this  area. 


{ 


WIND  SPEED  VARIABILITY  AND  STATISTICS 

Time  series  records  of  surface  wind  speed  and  atmospheric  pressure  are  displayed  in 
Figs.  1  through  4.  Each  plot  is  based  on  values  recorded  every  12  h,  at  noon  and  midnight, 
for  2  yr.  The  data  are  smoothed  with  a  uniform,  rectangular  window,  giving  a  3-day  (6-data- 
point)  average.  Visual  comparison  of  plots  with  and  without  smoothing  indicated  that  the 
smoothed  curves  are  better  able  to  reveal  the  variability  and  trends  over  the  time  periods  of 
interest.  Each  curve  shows  an  annual  seasonal  variability,  with  wind  speed  generally  greater 
in  winter  than  in  summer  and  atmospheric  pressures  generally  greater  in  summer.  It  is  also 
apparent  in  these  figures  that  variability  over  shorter  time  periods,  from  a  few  days  to  a 
month,  is  greater  than  the  seasonal  variability. 

Monthly  statistics  of  wind  speed  at  each  station  are  summarized  in  Table  1.  The  first 
column  for  each  station  gives  the  number  of  samples  available  each  month.  This  number 
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Fig.  3  —  Smoothed  time  series  record  of  wind  speed  et  Jen  Mayen, 
June  1974  through  Msy  1976 


Fig.  4  —  Smoothed  time  series  record  of  atmospheric  pressure  st  Jan  Mayen, 
June  1974  through  May  1976. 
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varies  partly  because  some  wind  speed  values  are  missing.  Values  of  standard  deviation, 
skewness  and  kurtosis  in  the  table  are  based  on  the  following  definitions  [1] : 


Standard  deviation  = 


N 


1/2 


Skewness 


2 No3 


2  /<*,  ~  X) 

Kurtosis  =  - 

2 No* 


3 

2 


(1) 


(2) 


(3) 


Monthly  values  of  standard  deviation  show  moderate  correlation  (correlation  of  0.70  at  Jan 
Mayen  and  0.65  at  Bear  Island)  with  monthly  average  wind  speeds,  a  tending  to  be  higher 
when  the  average  wind  speed  is  higher.  Values  of  skewness  and  kurtosis  show  relatively  large 
fluctuations  from  month  to  month. 


PROBABILITY  DENSITY  AND  CUMULATIVE  DISTRIBUTION  FUNCTIONS 

Probability  density  functions  and  cumulative  distributions  of  wind  speed  were  com¬ 
puted  for  each  station  and  are  shown  in  Figs.  5  through  8.  The  appearance  of  the  proba¬ 
bility  density  functions  is  similar  to  that  of  the  Rayleigh  distribution 


where  u  is  the  wind  speed  and  S  is  an  adjustable  parameter.  The  cumulative  distributions 
give  the  percentage  of  wind  speeds  less  than  the  abscissa  and  are  plotted  on  probability 
paper.  A  cumulative  Rayleigh  distribution,  given  by  1  —  exp(-u2/2S2),  was  fitted  to  the 
data  of  each  station,  and  the  corresponding  curves  are  indicated  in  Figs.  6  and  8.  The 
value  of  the  parameter  S  for  Station  1028  was  chosen  such  that  the  median  wind  speed 
is  the  same  for  the  fitted  curve  and  the  data.  This  value  is  S  *  5.65,  determined  from 
the  relation 


median  wind  speed 
(2  fin  2)1/2 


(5) 


5 


Fig.  5  —  Probability  denaity  function  of  wind  apeed  at  Bear  Ialand, 
June  1974  through  May  1976 


WIND  SPEED  (m/t) 


Fig.  6  —Cumulative  distribution  of  wind  speeds  at  Bear  Island,  June 
1974  through  May  1976.  Curve  corresponds  to  Rayleigh  distribution 
with  S  »  5.66. 


Fig.  8  —  Cumulative  distribution  of  wind  speeds  st  Jan  Mayen,  June 
1974  through  May  1976.  Cunres  correspond  to  Rayleigh  distribution 
with  S- 5.05  and  5  -  6.00. 
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The  resulting  computed  curve  agrees  well  with  the  data,  shown  by  the  symbols.  A  similar 
approach  for  Station  1001  gives  S  “  5.05;  the  resulting  computed  curve  shows  poorer  agree¬ 
ment  than  was  found  at  the  other  station.  The  cumulative  Rayleigh  distribution  for  S  =  6.00 
is  illustrated,  also,  for  the  purpose  of  comparison. 

Other  tests  that  can  be  used  to  measure  the  closeness  of  the  data  to  a  Rayleigh  distribu¬ 
tion  are  to  compare  values  of  the  coefficient  of  variation  (defined  as  the  standard  deviation 
divided  by  the  average),  the  skewness,  and  the  kurtosis.  Calculated  values  of  these  quantities 
for  a  Rayleigh  distribution  are  compared  to  the  measured  values  at  each  Station  in  Table  2. 

It  is  noted  that  the  calculated  values  for  a  Rayleigh  distribution  are  independent  of  the 
parameter  S. 

It  is  concluded  from  Table  2  as  well  as  from  Figs.  5  through  8  that  as  an  approximation 
the  distribution  of  wind  speeds  can  be  modeled  as  a  Rayleigh  distribution.  This  model  pro¬ 
duces  a  better  fit  to  the  data  at  Bear  Island  than  at  Jan  Mayen. 


FLUCTUATION  SPECTRA 

The  time  series  of  wind  speed  and  atmospheric  pressure  values  at  each  station  were 
analyzed  into  frequency  components  by  means  of  a  discrete  Fourier  transform.  The  resulting 
frequency  components  were  then  used  to  form  the  wind  speed  and  pressure  fluctuation 
spectra,  which  are  presented  in  Figs.  9  through  13.  Figures  11  through  13  are  based  on  a  2- 
yr  record,  data  taken  every  12  h,  for  a  total  of  1462  samples  of  each  variable.  Each  series 
gives  an  upper  spectral  frequency  of  1/24  cph.  Figure  10  is  based  on  a  6-month  record  from 
July  through  December  1974,  with  values  every  3  h.  The  higher  sample  rate  results  in  an 
upper  spectral  frequency  of  1/6  cph.  The  spectrum  in  Fig.  9  is  combined  from  both  the  2-yr 
and  the  6-month  records  and  has  the  broadest  spectral  range  available  from  the  data.  The 
two  spectral  portions  are  joined  at  a  frequency  of  approximately  0.026  cph,  the  lower- 
frequency  portion  coming  from  the  2-yr  record,  and  the  higher-frequency  portion  from  the 
6-month  record.  The  wind  fluctuation  spectra  in  Figs.  9, 10  and  12  have  been  smoothed 
with  a  sliding  1-3-&-3-1  frequency  window. 

The  largest  single  frequency  component  of  each  variable  is  the  annual  component.  For 
wind  speed,  the  peak  amplitude  annual  components  at  Bear  Island  and  Jan  Mayen  are 
1.3  m/s  and  1.7  m/s,  respectively.  (Note  that  the  ordinates  in  Figs.  9  through  13  give  the 
spectral  density,  which  is  the  mean  square  of  the  Fourier  frequency  component  divided  by 


Table  2.  —  Comparison  of  Measured  and  Rayleigh  Values 
of  Coefficient  of  Variation,  Skewness,  and  Kurtosis 


Rayleigh 

Sta.  1001 
Jan  Mayen 

Sta.  1028 
Bear  Island 

coeff.  of  var. 

0.52 

0.62 

0.49 

skewness 

0.32 

0.28 

0.31 

kurtosis 

0.12 

-0.01 

0.16 

10 
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the  fundamental  frequency.)  These  annual  components  are  relatively  small  as  compared  to 
the  total  variability  of  wind  speed  values.  The  variance  associated  with  the  annual  wind 
speed  component  is  0.88  (m/s)2  at  Bear  Island  and  1.44  (m/s)2  at  Jan  Mayen,  in  each  case 
much  less  than  the  corresponding  total  variances  of  14.1  (m/s)2  and  18.3  (m/s)2,  respec¬ 
tively.  Thus,  seasonal  variability  constitutes  a  relatively  small  portion  (approximately 
7  percent)  of  the  total  variance  of  wind  speed  at  each  station.  The  question  of  where  in 
the  spectrum  most  of  the  variability  is  found  is  considered  later.  Peak  amplitude  annual 
components  of  atmospheric  pressure  are  760  Pa  (7.6  mbar)  and  700  Pa  (7.0  mbar)  at  Bear 
Island  and  Jan  Mayen. 

At  high  frequencies,  the  wind  fluctuation  spectrum  decreases  at  a  rate  that  appears 
slightly  greater  than  f  ~ 1 ;  at  lower  frequencies  the  slope  is  much  less.  Transition  from  low  to 
high  fluctuation  frequency  for  the  wind  speed  appears  to  fall  in  the  spectral  band  near 
10"  2  cph,  corresponding  to  a  period  of  about  4  days.  The  6-month  spectrum  of  Fig.  10 
shows  a  moderate  peak  at  10"  2  cph. The  pressure  fluctuation  spectra  have  a  high  frequency 
slope  of  about  f~  3/2,  with  a  transition  to  the  high-frequency  asymptote  at  approximately 
0.008  cph,  corresponding  to  a  period  of  about  2  weeks. 

Figures  14  through  18  present  a  quantity  that  has  become  known  [2]  as  the  logarithmic 
spectrum,  each  corresponding  respectively  to  the  spectra  in  Figs.  9  through  13.  The  ordinate 
represents  the  product  of  the  fluctuation  spectrum  times  the  fluctuation  frequency,  plotted 
against  the  logarithm  of  frequency.  An  advantage  of  the  logarithmic  spectrum  is  that  the 
variance  for  any  portion  of  the  spectrum  is  proportional  to  the  area  undo:  the  curve,  while 
preserving  the  advantages  of  using  a  logarithmic  frequency  axis.  Use  of  the  logarithmic 
spectrum  aids  in  the  recognition  of  major  spectral  contributions  in  proportional  frequency 
bands,  rather  than  in  bands  of  fixed  bandwidth.  The  logarithmic  spectra  in  Figs.  14  through 
18  are  computed  from  fluctuation  spectra  that  were  first  smoothed  by  a  triangular  sliding 
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Fig.  18  —  Logarithmic  pressure  spectrum  at  Jan  Mayen,  June  1974  through  May  1976 


frequency  window,  averaging  over  five  successive  values  with  relative  weightings  in  the 
ratio  1-3-5-3-1.  Figure  14  presents  combined  results  based  on  the  2-yr  as  well  as  the  6-month 
data  records  in  the  same  manner  as  Fig.  9.  The  logarithmic  wind  spectra  show  that  most  of 
the  variability  exists  in  a  broad  spectral  region  whose  roughly  defined  limits  correspond  to 
periods  of  from  2  to  10  days.  Figure  15  shows  a  peak  at  0.01  cph,  based  on  the  6-month 
record,  but  this  peak  is  not  apparent  in  the  2-yr  records.  The  pressure  records  show,  similarly, 
that  most  of  the  spectral  energy  exists  in  a  broad  spectral  band  whose  roughly  defined  limits 
correspond  to  periods  of  from  5  to  20  days. 

In  performing  the  Fourier  analysis  of  the  2-yr  time  series  of  values  of  wind  speed  or 
pressure,  it  was  regarded  as  necessary  to  use  the  entire  2-yr  sequence  of  1462  values,  in 
order  that  the  annual  variation  coincide  with  a  harmonic  (the  second)  of  the  fundamental 
analysis  frequency.  Thus,  use  of  the  available,  conventional  FFT  algorithm  was  not  possible 
because,  based  on  a  time  length  other  than  the  full  2-yr,  it  would  tend  to  smear  out  the 
annual  spectral  contribution.  In  order  to  replace  the  FFT  while  keeping  computation  time 
small,  a  separate  discrete  Fourier  transform  subroutine  was  written,  based  on  a  conventional 
approach  but  exploiting  the  high  speed  vectorizing  capability  of  the  TI-ASC  computer.  The 
subroutine  can  be  adapted  to  any  number  of  data  points  and  is  presented  in  the  Appendix. 
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DISCUSSION  AND  SUMMARY 

The  result*  of  several  early  studies  of  wind  speed  fluctuations  are  summarized  by 
Van  der  Hoven  [3]  in  a  report  that  covers  an  extraordinarily  broad  spectral  range.  His  aim 
was  to  cover  a  sufficiently  broad  part  of  the  spectrum  "so  as  to  include  all  important  con¬ 
tributions  to  the  total  variance.”  Van  der  Hoven  identifies  two  major  spectral  peaks,  one 
centered  at  a  period  of  4  days  (with  a  spread  from  about  2  to  7  days)  attributed  to  the 
passage  of  synoptic  scale  weather  patterns,  the  other  at  a  period  of  about  1  min  and  related 
to  the  micrometeorological  turbulence.  The  reported  peak  at  4  days  has  a  logarithmic  spec¬ 
tral  level  of  about  5  (m/s)2.  His  results  are  generally  consistent  with  the  results  presented 
here. 


The  existence  of  a  4-day  periodicity  has  been  reported  also  in  regard  to  oceanographic 
variables  whose  response  appears  forced  by  atmospheric  variation.  Thus,  Groves  [4]  observed 
a  4-day  spectral  peak  in  time  series  of  records  of  sea  level  in  equatorial  regions  of  the  Pacific 
Ocean,  correlated  with  similar  fluctuation  in  wind  speed  data.  Christensen  and  Rodriguez 
[5]  note  a  5  to  7-day  oscillation  in  ocean  currents,  supposedly  driven  by  atmospheric  varia¬ 
tion,  off  the  coast  of  Baja,  California.  Similarly,  Wunsch  and  Wimbush  [6] ,  and  Diiing, 
Mooers,  and  Lee  [7]  identify  a  broad  spectral  peak  with  periods  from  4  to  7  days  in  current 
fluctuations  off  the  Florida  coast,  correlated  with  atmospheric  variability. 

Two  principal  conclusions  are  drawn  from  the  present  analysis  of  wind  speed  data  from 
the  two  stations  in  the  Norwegian  and  Barents  Seas,  as  well  as  from  earlier  studies  of  wind 
speed  fluctuation. 

1.  Wind  speed  values  considered  over  a  sufficiently  long  time  period  follow  a  proba¬ 
bility  distribution  that  can  be  modeled,  at  least  to  a  first  level  of  approximation,  as  a  Ray¬ 
leigh  distribution. 

2.  Wind  speed  fluctuations  contain  strong  components  at  frequencies  of  one  cycle  per 
year  and  its  harmonics.  Most  of  the  fluctuation  power,  however,  occurs  in  a  fairly  broad 
spectral  band  with  periods  between  2  and  10  days.  Because  of  this  large  spectral  concentra¬ 
tion  of  fluctuation  power  at  from  2  to  10-day  periods,  it  follows  that  ocean  acoustic  experi¬ 
ments  whose  purpose  is  to  characterize  seasonal  variability  of  weather  related  phenomena, 
such  as  ambient  noise,  should  be  conducted  over  time  periods  sufficiently  long  to  allow 
averaging  over  these  synoptic  variations. 
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Appendix 

FOURIER  ANALYSIS  OF  WIND  SPEED  DATA 

The  wind  speed  can  be  expressed  as  a  Fourier  Series 

N  N 

Wk=Ao+^L  A"  COS  2  ™  f  +  S  Bn  sin  2  ™  f  • 
n=*  1  °  n-1  ° 

where, 

Ta  is  the  basic  period  of  the  data,  and 
N  is  the  highest  term  in  the  Fourier  Series. 

In  this  case  the  data  were  taken  at  12-h  intervals  over  a  period  of  2  yr  giving  Ta  =  2  yr  = 
1462  X  12  h  and  N  =  1462/2  -  731. 

Furthermore,  time  t  changes  in  discrete  steps,  so  t  can  be  replaced  by 

t  =  fet0,  (t„  =  12h) 


giving 


N  N 

Wk  =  Ao  +  ]C  A »  C°S  2  nnk/U 62  +  S  Bn  Sin  2  ™fc/1462- 
n»l  n=l 


By  multiplying  both  sides  of  this  equation  by  either  sin  (2  Trnfe/1462)  or  cos  (2  nnk/1462), 
summing  over  fe  from  1  to  1462,  and  using  orthogonality  rules,  one  finds 


A 


n 


1462 


*=1 


2  nnk 

W.  cos - —  ,  and 

*  1462 


B 


n 


1 

731 


1462 

E 


*=i 


Wk  sin 


2  nnk 
1462  ’ 


A0  is,  of  course,  the  average  of  W* . 


The  calling  sequence  for  the  Fourier  transform  subroutine  is 

CALL  ANBN 
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where 

IV  =  a  Real  *4  array  containing  the  time  series  Vfk> 

K  =  an  Integer  *4  variable  defining  the  number  of  data  points  in  the  time  series  Wk, 

A  =  a  Real  *4  array  containing  the  An  output  elements  of  the  transformation.  The 
length  of  the  array  is  N, 

B  =  A  Real  *4  array  containing  the  B„  output  elements  of  the  transformation.  The 
length  of  the  array  is  N,  _ 

C  =  a  Real  *4  array  containing  the  \J A^  +  B%  for  each  n .  The  length  of  the  array 
is  N,  and 

N  =  an  Integer  *4  variable  defining  the  total  number  of  terms  in  the  Fourier  Series. 

The  total  length  of  the  subroutine  is  00009300.  It  is  compiled  at  the  K  level  of  the  NX 
compiler.  The  program  calculates  values  of  the  sin  and  cos  using  their  recursive  feature  to 
compute  the  N  terms  in  the  Fourier  Series.  This  speeds  up  the  calculation.  For  this  case 
where  K  =  1462  and  N  =  731  the  total  central  processor  time  was  3.996  s.  The  program 
uses  the  following  external  routines:  COS,  SIN,  DOTPRD,  SQRT  and  FLOAT. 
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ooci 

SUBROUTINE  ANBN<W,K, A,B,C,N) 

99C2 

DIMENSION  W(K/200C/)»«(N/lCOO/)»B(*/1000/)fC(N/1000S)« 
1VCSSC2C0C )fVSNN(2000)«VC(2000)»VS( 2000) 

0<*C3 

PHASE*<2.*3.14159265)/FL0AT(K> 

09C4 

CON-1.0/FLOAT(N> 

0915 

00  5  1=1, K 

00C6 

VCSSCL)*COSCPHASE*L) 

00c  T 

5 

VSNNC L)ZS INC PHASE* L) 

0938 

DO  7  1=1, N 

0909 

18 EGZI 

0910 

INC*  I 

0011 

J«0 

0)12 

10 

CONTINUE 

0913 

00  9  M=IBEG, K , INC 

0014 

JZJ*1 

0915 

NTOP*N 

0916 

8CCJ)*VCSSCN3 

0017 

9 

VSC J)*VSNN(N) 

0018 

IFCJ.EQ.K)  GO  TO  11 

0019 

IBEGzCNTOP*INC)-K 

0920 

GO  TO  10 

0921 

11 

AC1)ZCON*OOTPROCU,VC) 

0022 

8CI)=C0N*D0TPR0C«»VS) 

0923 

CSQ=ACI)*ACI>+BCI)*BCI) 

0)24 

C<n*SQRTCCSQ) 

0925 

7 

CONTINUE 

0026 

RETURN 

0027 

ENO 

